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Abstract 

For the first order transition of the Ising model below T c , 
Isakov has proven that the free energy possesses an essential sin- 
gularity in the applied field. Such a singularity in the control 
parameter, anticipated by condensation theory is believed to be 
a generic feature of first order transitions, but too weak to be 
observable. We study these issues for the temperature driven 
transition of the q states 2D Potts model at q > q c = 4. Adapt- 
ing the droplet model to this case, we relate its parameters to the 
critical properties at q c and confront the free energy to the many 
informations brought by previous works. The essential singular- 
ity predicted at the transition temperature leads to observable 
effects in numerical data. On a finite lattice, a metastability do- 
main of temperatures is identified, which shrinks to zero in the 
thermodynamical limit. 
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1 Introduction 



The transition between the vapour and liquid phases of a material has been 
for a long time the prototype for first order transitions. Below a critical 
temperature T c , the density, that is the derivative of the pressure V with 
respect to the chemical potential /i, is discontinuous at \i = fj, t , where the two 
phases coexist. The vapour phase is stable below /x t , the liquid phase above. 
As T — * T c , the transition becomes continuous, the density fluctuations 
diverge, leading to the well known phenomenon of critical opalescence. Van 
der Waals theory of condensation predicts these equilibrium properties, and 
in addition lead to the existence of metastable states, vapour above /it and 
liquid below, experimentally accessible during finite times. In this theory, 
they correspond to the possibility of analytically continuing the low and 
high density branches of the isotherm function V(fi) \t beyond /it ■ 

All of this can be found in standard textbooks, together with other 
methods leading to similar descriptions, as well as applications of the same 
ideas to many systems of physical and/or theoretical interest, which undergo 
first order transitions. One of the simplest and most extensively studied 
system is the Ising model. There instead of the pressure, one considers 
J-{h) \t, the free energy as a function of the applied magnetic field h at 
fixed T. Below T c , the field driven first order transition is located at ht = 0, 
and manifested by a discontinuity of the magnetization. As T — > T c , the 
magnetization vanishes, and its fluctuations measured by the susceptibility, 
diverge with a critical exponent which is known in 2D. In this particular case 
however, it has been mathematically proven by Isakov Q that Van der Waals 
theory cannot be true: each of the ordered phase free energies possesses 
an essential singularity at h = ht, which makes its analytical continuation 
ambiguous in the metastability domains where the magnetization remains 
opposite to the applied field. 

In fact, such a singularity in the control parameter at the transition 
was anticipated a long time ago from the droplet theory of condensation in 
the version developped by Fisher H, who attributes to Mayer ||] the first 
suggestion in this direction. Shortly after, Langer H applied the droplet 
model to the Ising case, emphasizing that the expected singularity was so 
weak that it could not be experimentally detected. The reader will find in 

||] many references to earlier and related work, and inspiring discussions 
on the physics behind. Since then, the case has been reexamined many 
times, with particular concerns about the definition and properties of the 
so-called droplets, or clusters. Steps in this line of research can be followed 
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in, e.g., [|5|, ]6|, [jj. Its intimate connection with metastability and dynamics 
at first order transitions has been for a long time a matter of investigations 
and debate f§ §, and still is an active field 0, [□], 0, 0, 0, |l5|. 

The basic idea of the approach is that the fluctuations of the system in 
a stable phase are due to the appearance of droplets of one or several other 
phases, whose size and number are controlled by the value of the driving 
parameter. The statistics of these droplets is described as a Mayer's cluster 
expansion of the thermodynamical function at hand, and the analyticity 
of the latter follows from the convergence properties of this expansion. A 
transition occurs at a point of parameter space where the energy/entropy 
balance favours large clusters of the wrong phase, and the cluster expansion 
is shown to diverge at such a point, which thus is singular. 

In this paper, we investigate similar issues within the bidimensional, te- 
states, Potts model [16] for which we adapt the droplet picture, much in the 
spirit of [^], using the fact that the transition, here driven by temperature 
and being first order for q > 4, becomes continuous as q — > q c = 4 where the 
latent heat and the interfacial tension vanish. The Potts model is interesting 
per se, as illustrated by the many studies which have been and still are 
devoted to it, and it seems thus worthwhile to learn more about it. As 
we shall see, the droplet point of view allows one to put together in a very 
compact way many of the results accumulated over years, either analytically 
or by high precision simulations [17, 15] at the transition. It is the case for 
the relationship between properties of the free energy along the first order 
transition line q > 4 and at its critical end point q c , a relationship which 
was evidenced in a series of papers by Bhattacharya et al. |l^, [2(], pj]] . 
This connection was recently studied in a completely different framework by 
Cardy [22] who relates the emergence of a critical q c to branching properties 
of the Potts interfaces between different ordered phases. A better knowledge 
of the analytic structure of the free energy close to a first order point may also 
help in other domains of current interest, such as the dynamics, as already 
mentioned, and the role of disorder for the phase diagram of statistical 
systems. 

Duality relates the disordered phase, stable at high temperature, to the q 
ordered phases, stable at low temperature. We will mainly focuss on the dis- 
ordered phase. The paper is organized as follows. In Section (2), we adapt 
the droplet formulation to the Potts case, and constrain its parameters by 
requiring consistency, for q values above 4, with the known critical properties 
at q c . In Section (3), we first show that, in a pure phase, the order n energy 
cumulants f n of the model compare successfully with those measured in [[D]] 
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for n < 8 and q = 10, 15, 20. Only one, g-independent, parameter is needed 
to reproduce all ratios f n IH given fi at each q. This scaling property allows 
one to construct a universal (i.e. (/-independent) function 4>(z), which repre- 
sents the pure phase free energy suitably rescaled, z measuring the distance 
to the transition point in units of the inverse temperature /3, also rescaled. 
The function <fi has a simple integral representation from which we derive 
its analytic properties. It has an essential branch point singularity at z = 0. 
We then turn to applications. The associated internal energy distribution 
for a finite square lattice is studied in Section (4). In particular, rescaling 
the internal energy and the lattice size leads to a universal distribution from 
which any g-dependence has disappeared. We emphasize consequences for 
the distribution of the singular structure of the free energy at the transition. 
Specific finite size effects are predicted, although by assumption the ther- 
modynamical limit of the free energy density has been taken, and compared 
with those actually observed in a numerical simulation performed at q=9. In 
Section (5), we define the disordered distribution below the transition tem- 
perature by reweighting that defined above. From its finite size properties, 
we determine a size dependent spinodal value (3* such that the system is 
metastable in the interval (3* > (3 > (3f Inside this interval, which shrinks 
to zero in the thermodynamical limit, we define a free energy for metastable 
states . A summary and conclusions are proposed in a last section. 

2 The Potts model and its droplet description 

We first recall a few basic definitions and properties relevant for our dis- 
cussion of thermodynamical quantities associated with the ferromagnetic 
g-states Potts model. We next summarize the standard droplet descrip- 
tion of condensation, and finally adapt it so as to incorporate the known 
properties of the model at q > 4 . 



where < ij > denotes the pairs of nearest neighbours on a square lattice of 
area A = L 2 and o"j one of the q possible values of the spin variable at site 



2.1 The Potts Model 



The model is defined JT^j through the Hamiltonian 





<ij> 
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i. In the rest of the paper, the energy per link J is set equal to 1. 
The partition function of the system is 

Z A {(3) = 5>xp(-/3tf(a)), (2) 
W 

= J2n A (E)exp(-{3AE) (3) 

E 

where (3 is the inverse temperature, E G [—2, 0] the energy density per site, 
and £1 A the corresponding number of states. It is equivalent to say that the 
energy probability density is 

PpAE) = V A (E) exp(-/3AE). (4) 

For convenience, we shall call "free energy" the quantity 

F A = jlog(Z A ), (5) 

while strictly speaking F A is —1//3 times the standard free energy density. 
Its thermodynamical limit F is its limit as A — ► oo. The finite size micro- 
canonical entropy density S A (E) and its thermodynamical limit S(E) are 
defined by 

Sa(E) = ^log(OA), (6) 

= F A (J3) + jlog{P PtA (E)) +(3E, (7) 
S(E) = UmS A (E), A^oo. (8) 

A numerical simulation provides the energy distribution Q A e~^ AE up to a 
numerical factor, and thus the entropy S A up to an additive constant. 
A review of many of the known properties of the Potts model can be found 
in 1 22]. A transition occurs at an inverse temperature (3t given by p4j] 

Pt = log(^ + l). (9) 

It is first order for q > 4, and second order at and below the end point 
q c = 4, where it is characterized, e.g., by the critical indices a and v for the 
specific heat and the correlation length respectively . For future reference, 
we recall that 

a = 2/3, 

v = 2/3. (10) 
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In the thermodynamical limit, the system is either in the disordered phase, 
for P < fa, or in one of the q degenerate ordered phases for P > p t . When 
necessary, a quantity referring to an ordered or to the disordered phase will 
receive a superscript o or d respectively. So (resp. F^) denotes the 
ordered (resp. disordered) free energy, a well defined function of P for (3 > (3 t 
(resp. < Pt). The transition point is the P value where the two free energies 
are equal: F^°\p t ) = F^ d \p t ). For q > 4 where the transition is first order, 
one expects on very general grounds that at Pt, F^ (resp. F^) has finite 
left (resp. right) derivatives to any order n with respect to p. They give 
(— l) n times the cumulants f°> d of the internal energy associated with each 
of the pure phases o or d. We shall be interested in the analytic continuation 
of the free energies in the complex P plane at q > 4 fixed, which will consist 
in giving a meaning to the formal expansion: 

oo 

F M (P) = ]T {-l) n f° n d {P - PtT/n\ . (11) 

n=0 



The two functions F are related to each other by duality [23]. If 

(exp(^) - 1) (exp(/3) - 1) = q, (12) 

then 

F( )(^) = F(°)(/3)-2 1og((exp(/3)-l)/^) . (13) 

We choose to study F^ for concretness, and consider the d —* o transition 
as a condensation process where droplets of aligned spins inside a disordered 
bulk tend to grow as P approaches Pt from below. The droplet model consists 
in describing the system at equilibrium as a statistical distribution of non 
interacting droplets. 



2.2 The Droplet Picture 

We closely follow Fisher The transposition from vapour condensation 
to the d — > o Potts transition at q > q c produces the disordered phase free 
energy F^ d \p) as a function of the inverse temperature. A droplet of size 
£ (its area) is a connected domain of £ sites with the same spin value. The 
model assumes that it has an effective perimeter scaling as £ a , where the 
exponent a, smaller than 1 (no fully ramified clusters), but possibly larger 
than the geometrical value 1/2 [pq] , accounts for the appearance at fixed £ 
of many different shapes and topologies. This is controversial and calls for 
comments to be made later. The free energy can then be written: 
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F^(fi) = cJ2r T x f °y e (14) 
e=i 

y = exp(/?-A) (15) 

x = exp(— (J) (16) 

with 

u > 0, 1/2 < a < 1, r > 0. (17) 
The parameter c fixes the normalisation and ^ r is a correction to the area 



and perimeter dependences EL 26]. A priori, the parameters u,a,r and c 



are functions of q, a question which we examine in the next subsection. Eqs. 



(14-16) are interpreted as follows. The free energy is the sum over £ of 
contributions coming from all clusters of size I. The quantity y plays the 
role of the activity in gas condensation or of exp (—h) in the Ising problem, 
and ijj denote an energy density per unit of effective perimeter (effective 
interface tension between the interior and the exterior of a droplet). 

Each term in the above sum is proportional to the probability of an or- 
dered cluster of size I in the system. Due to inequalities (|l7|), the factor yr 
implies that for (3 < fit the probability of arbitrary large clusters is expo- 
nentially small. The disordered phase is stable there, the series and all its 
derivatives with respect to fi, converge. The converse is true for fi > fit- 
It follows that the disordered free energy is analytic in the half-plane 
Re(fi) < fit and has an essential singularity at fi t . 



As a side remark, we note that in its Kasteleyn-Fortuin formulation [27], 
the Potts model is a model of satisfied or unsatisfied links. The temper- 
ature dependence of the partition function is carried by weigths y l , with 
y = (exp(/3) — l)/y/q, t being the number of satisfied links forming a con- 
nected cluster. The droplet picture applied to the link representation leads 



to equations similar to (14,15) with y replaced by y. Close to the transition, 
y ~ y ~ 1 up to order (fi — fi t ) terms, so that using y or y makes no differ- 
ence for the leading singularity: It has the same structure in b = logy and 
in/?- fit. 

Following Fisher Q, we now consider the system close to the critical 
point q c , and fix the parameters u, a and r. 
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2.3 Fixing the Parameters from the Critical Point 

We go on with the disordered phase free energy, and from now on omit the 
d superscript. The energy cumulant of order n is given by 

fn = (-l)"|^l*=l ( 18 ) 



with F given by Eq.(|l4|). At least for large enough orders, the series obtained 
for f n can be given a closed form, by replacing the discrete sum over I by 
an integral M, m. The result is 



/„ = c (_i)« / d£exp(-u;t)r- T (19) 
Jo 

c(-ir 



{n-T+l)/a 



r((n-r+l)/o-). (20) 



Hence all the derivatives of F exist at the transition, but due to a < 1 the 
convergence radius of its Taylor series (|Tl| ) in (j3 — /3t) is zero, a characteristic 
situation for an essential singularity at 0f It also follows from (^) that the 
scale in f3 — fit relevant for the behaviour of F close to the transition is 

Po = ^. (21) 

On another hand, from an expansion to order 10 in l/-y/q of the free energy, 
it was empirically found in 2C, 21], and confirmed by expansions recently 



pushed to order 23 [28], that at least in a "low" q domain, in fact extending 
up to q ~ 30 or more, one has approximately 

h « e, (22) 

hlh oc e'\ (23) 
where £, the correlation length of the Potts model, is exactly known |2^. R^] 



and grows extremely rapidly as q decreases towards q c . In (l^, 20], a tentative 



explanation for such behaviours was that at q c , £ and the most singular part 
Fsing. of F are respectively proportional to \(3 — (3t\~ u and |/3— (3t\ 2 ~ a , so that 



Eqs.(22,E3) hold there, due to fllOD, and thus perhaps in some neighborhood 



of the critical point. Together with fl20j, pi[), this invites us to postulate 

u = /// 2 cxl/£, (24) 
a = 2/3, (25) 
r = 7/3, (26) 
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with / approximately constant with respect to q. 

Note that if u in (14,^) is interpreted as an effective interfacial tension 
between the interior and the exterior of a droplet, its postulated behaviour 
(|24|) fits nicely that of the true o—d interface tension X, known to be equal to 
l/2£ p0[ , and also determined numerically via a study of phase coexistence 
at (it [HI]- We remark that at q = q c the same behaviour X£ — > est. as 
/? — * P c follows from the scaling relation fj, = v, fj, being the Widom exponent 
for the interface tension [|31{|. Why it also holds for cut; remains however 
unclear since whereas X is an energy density per unit length along a straight 
boundary, u> is here related to an effective boundary whose size scales as a 
length to the power 2a = 4/3 (|25|), 

As discussed in details by Binder [|[] for the Ising case, taking a ^ 1/2 
is quite questionable, and in fact incorrect if one insists that the droplets 
of the model are the standard geometric clusters which can be drawn on a 
lattice. This is confirmed by the exact results of || about the probability 
of such clusters, and of Isakov ||] on the cumulants of the Ising model 
at sufficiently low temperature. Because their derivation is essentially of 
geometrical nature, we see no reason why they should not apply to the Potts 
case. A physical interpretation of a > 1/2 has been proposed by Binder. For 
the time being, we take a = 2/3 as an effective parameter describing some 
neihborhood of the critical end point, and show that anyway it accounts for 
many observed facts. This point will be discussed again in our conclusions 
(Section (6)). 



3 The Cumulants confronted to Numerical Data. 
The Free Energy 



The arguments borrowed from jjl], 20 1 to guess fs/f2 Eq.(23) were fur- 
thermore applied in these references to cumulants of any order by successive 
derivations of F s i ng . cx {(3t—(3) 2 ~ a and the subsequent replacement of (/3j— (3) 
by a constant times As a result, the series (11) converged, providing 

an ansatz for |2l]] analytic up to a value j3* > fit, which could be inter- 
preted as the end of a metastability region (spinode) , in clear contradiction 
with any droplet approach. 
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3.1 The Potts model cumulants 



Let us compare the model predictions for the energy cumulants with the 
results of the high accuracy simulation of Janke and Kappler [17], who give 
the pure phase f n up to order 8 at Pf We are concerned with n > 2 only, i.e. 
with those cumulants which diverge at the critical point. Anyway, /o and f\ 
are exactly known p4|], and irrelevant for the study of a single phase. Since 
the normalisation constant c remains undetermined and since f% is supposed 
to carry most of the ^-dependence of u> (|24|), we rewrite fl20|) as 



fn 
h 



(-iy 



h\ {n - 2),a T((n 



r + l)/a) 



r((3-r)/a) 



n > 2, 



that is, if the conjectures (25, 26) are used 



fn 
T2 



(-l) n 5 n r(3n/2-2), 

/ 2 \(3«/2-3) 
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(27) 

(28) 
(29) 



This is the basic equation of our paper. It determines all the / n 's, n > 2 as 
a function of fi at the price of one free parameter / only. We fitted Eq. ( p8| ) 
to the data [17], fixing fi at the central value measured for each (?=10, 15 
and 20, n varying from 3 to 8. We obtained a x 2 /d.o.f. equal to 1.3 for 17 
degrees of freedom with 



/ = 0.295 ± 0.003 



(30) 



The excellent agreement is illustrated in Fig.[l] where the numerical data 
JL7[ for (— l) n / n /(/2 S n ), n < 8 are compared to T(3n/2 — 2) (continuous 
line) and plotted against the order n. Also indicated are the orders of 
magnitude estimated for n = 9 and 10. In order to appreciate whether the 
values ((2^, HH) of a and r are actually requested by the numerical data, we 
have repeated the fit letting also these parameters free. The new fit requires 



a 

T 
f 



0.67 ±0.01, 
2.34 ±0.10, 
0.28 ±0.02, 



(31) 
(32) 
(33) 



with x 2 /d.o.f. = 1 for 15 d.o.f.., in perfect agreement with our conjecture. 
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This nicely confirms the existence of a large domain of "low" q values 
where the fluctuations of the system are strongly influenced by the critical 
properties at q c . The Fisher droplet model, supplemented by our conjecture 
that the whole g-dependence is embedded in /2, gives them an economical 
and accurate description. We note that the ansatz |HJ] predicts much smaller 
values for the highest order cumulants, missing /§ by a factor about thirty. 
These results not only support the existence of an essential singularity in 
the free energy, but also show that this singularity does lead to detectable 
effects. We now proceed to construct the free energy explicitly . 



3.2 The Disordered Phase Free Energy 

The free energy can be reconstructed either directly from the droplet for- 
mulation (14), replacing the sum by an integral as done in Q, or from the 
asymptotic series (|ll]) with f n given by ( |27|) or (|28|) and the T-function re- 
placed by its integral representation (Borel resummation) . In view of the 
results of the previous subsection, it is natural to introduce the following 
rescaled free energy cfi, expressed as a function of a rescaled temperature z, 

Hz) = (^) 2/CT j 2 (F((3)-F(p t )-f 1 ((3-f3 t )), (34) 

z = -(/3-A)Cf 2 //) 1/<J . (35) 

With this choice of normalizations, we have <p(0) = (j>'(0) = and <ft"(0) = 1. 
For generic values of a and r, <t> reads 

m = W - r)M Jo 
and specializing to the values ( ^5[ |2^ ) of a and r, we finally get 

*(*) = f^e-^e-^-l + zt^). (37) 

This function is obviously holomorphic for Re z > 0, which includes the sta- 
bility region @ < of the disordered phase. It can be analytically continued 
in the complex z-plane by deforming the integration contour in t. The only 
singularity is a branch point at z = 0, easy to characterize. The contribution 
to <j) from any finite t interval [0, to] is an entire function, as well as that 
from to to oo of the linear part (in z) of the integrand. Hence the singular 
part of 4> is that of 

°°^e-^ 3/2 , (38) 
to t 



f ii;T | IW e-(e-"--l + ^), (30) 
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for any positive to- Next, as z is continued to z = \z\exp(±iTr), one may con- 
tinuously move the contour [to, +00] to [fo,T*o°] leading to a discontinuity 
along the negative real axis 

A = 



In the following, we define 4>{z) by ([37]) in the whole complex plane cut along 
the negative real axis, on which <p acquires an imaginary part. For \z\ small, 
the above integral can be estimated by steepest descent, giving 




where the function I(z), easy to evaluate numerically, goes smoothly to one 
as z goes to along the negative real axis. All this is quite similar to the 
analysis performed by Langer M for the Ising case, where the magnetic field 
plays the role of z. Whenever a precise evaluation of <j>(z) is needed, one can 
either numerically integrate the analytic continuation of ( |37| ) using suitable 
contours, or use its large \z\ expansion, which is found to be: 

<f,(z) = 7/6 - 3/4 -log(z)/3- 2^ + 5 (42) 
S = \z^ £ ^-^(ti), (43) 

where 7 is the Euler's constant. The derivation of this expansion can be 
obtained by a Sommerfeld- Watson method. One starts from the formal 
series expansion of <j) in z n , truncated at some large n value N. The sum 
J2o ( — l) n fln is replaced by the integral l/(2i) § dv a v j sin(7r v) over a contour 
encircling all integers from to N. After sending TV to 00 the contour is 
deformed so as to encircle the poles of the integrand situated on the negative 
real axis, and the residue formula is finally applied. A similar expansion for 
the second derivative of <j) with respect to z is much easier to derive starting 
from 

= r die-*'** 3 ' 2 , (44) 
J 

= z- 2 / 3 dtexp(-t/z 2 / 3 -t 3/2 ), (45) 
Jo 



\z\+ie)-<f)(-\z\-ie), (39) 
^e-*- rf3/ . (40) 



fo 
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and expanding the latter integral in powers of z~ 2 ^ 3 which yields a conver- 
gent series in this variable [32|. 

As we have seen, the closed form ( |37| ) of the thermodynamical free energy 
incorporates the information acquired upon the ~ 10 first energy cumulants 
in simulations performed on large but finite lattices during finite (Monte- 
Carlo) times. That it makes sense is justified by exact results [34], which, in 
short, state that on a finite lattice, 1/A log(Z^) gives the thermodynamical 
limit up to corrections which are exponentially small in the linear size 
of the lattice. In contrast, the next section will illustrate that the energy 
distribution Pa(E) on a finite lattice exhibits strong finite size effects, spe- 
cific to the singularity structure of the thermodynamical free energy, and 
corroborated by a numerical simulation at q = 9. The procedure followed 
to study Pa is similar to that used in [21]. 



4 Finite Size Effects and Scaling Properties in the 
Internal Energy Distribution 

4.1 The Probability Density for the Energy 

Given the rescaled free energy 4>(z), the probability density for observing 
a lattice averaged energy density E in a supposedly disordered bulk phase 
at /? = /?£, on a lattice of area A, is obtained by inverse Laplace transform 
using Eqs. (§, |, |). We find 

A rz+ioo 

PfrAt) = o± / dzexp(A r {<j>{z) -ez)), (46) 

Z17T J z—ioo 

where the rescaled energy, area and inverse temperature e, A r and z respec- 
tively are defined as 

1 / f \ 1 / 2 

" = /(/) ( E - E{d) ^ ^ 
A r = f{£) A i (48) 

Z = '{jT 2 ^'^- (49) 

The symbol (or E^) represents the exactly known internal energy of 
the disordered (or ordered) phase in the thermodynamical limit p3[. In 
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the rest of this paper, and depending on the context, we will use either 
the physical E, A and (3 variables, or their rescaled forms e, A r and z. For 
economy of notations however, we will often keep the same name for func- 
tions of them. The variable E (or e) is considered as continuous, (3 (or z) 
as unlimited in the imaginary direction. Eq.(|47|) maps the physical region 



E G [—2, 0] over a large interval in e, the more so f 2 is large (see Eq. fl47|), 
that is q close to 4. Since on large lattices P is very sharply peaked around 
e = 0, it will be justified (and convenient) to consider also e as unlimited, 
unless specified. Due to the analyticity properties of 4>{z), the integral (|46|) 
does not depend on z > 0. 

With the above prescriptions, one can verify that P^ tjJ 4(e) is a proba- 
bility density, normalized to one with respect to integration over e, that 
< e >= 0, < e 2 >= 1, and more generally that the cumulants of e for the 



density Pp tt A have exactly the values assigned by Eqs.(|34l, 35, £^), a very 
powerful check for numerical integration of (|4^). Any explicit ^-dependence 
has disappeared, and the effective area A r is the only external parameter. 
Its relation (^) to the physical area A predicts the scaling property that 
different q- values lead to the same Pp u A{t) if they are measured on lattices 
of linear size proportional to fa (i.e. approximately to £ according to (p 



which sounds reasonable). We note that the largest L values used in [17] for 
q = [10, 15, 20] were precisely chosen roughly proportional to £. 

This scaling behaviour is independent of the particular form of (f)(z), 
and we turn to the study of those properties of the distribution which are 
specific to our construction. The results of this section are all summarized 
through their application to the case q = 9 with L =80, 120 and 200, shown 



in Fig. [2]. The rescaling of A = L to A r (48) requires the knowledge of /2 



at q = 9. There the correlation length is 14.9 [pOl, and we evaluate f2 from 



its value at q = 10 by (§2|), finding f 2 (q = 9) = 12.8. The quantity 
shown in Fig. [2] is 

qPt M) = ^o g (^#) . (50) 



Ar V VA- 



The continuous lines correspond to the numerical integration of (46), and 



the other ones to various analytical approximations described below. 
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4.2 The Distribution above the Disordered Peak, e > 0. 

For e > 0, the saddle point equation for the integral (|46|), namely 

d4 _ 

dz 



(51) 



has a unique, non negative solution z s (e). Taking z = z s in the definition 
(f46|) makes the numerical integration easy. The saddle point estimate of 
Pp t ,Ai valid at large A r , is 



P 



Pt,A 



! A r <p"( z) 
2n 



exp 



d<j> 
dz 



z s {e), e > 0. (52) 



In this approximation, the thermodynamical limit of qp t A (50) exists for 
any non negative e, and defines, up to a linear term in e, a non trivial, 
concave, microcanonical entropy density S(E). Moreover, for e large, an 
estimate z as of z s can be obtained, by solving Eq.(|5l]) analytically when the 
two leading terms only of the large \z\ expansion ( |42"| ) are kept. Then 4>(z as ) 
and (j) {z as ) are computed at the same order, and plugged into fl5^), avoiding 
any integration. The corresponding algebraic expressions are quite lengthy, 
and not here. They are used for Fig. [2], showing that the approximation 
is very good in practice shortly after e=l. Farther away, these expressions 
show that <Z/3 ti oo( e ) become proportional to — (e + 2y / 7r) 4 . This means an 
extremely fast decrease of P(e) on the right hand side of the peak, much 
faster than any Gaussian. 



4.3 The Distribution below the Disordered Peak. Scaling at 

e < 

The side e < is more unusual, its properties reflecting the tendancy of 
the disordered system to order. As e passes through 0, the solution of 
( |5l|) disappears from the first Rieman sheet, and we need to optimize the 
integration contour in ( |46| ) differently. Given e, we move it to a contour 
z = x ±iy(x,e) passing through and determined by the condition 

Im((f)- ez) = 0. (53) 

In practice, because Im(4>(x)) is never large before the integrand in(|46|) gets 
exponentially damped at large \x\, we expand <f>(x + iy) in y to second order 
(note that eft is analytic and thus expandable at any non zero x), and solve 
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( p3D for y as a function of 4>(x) and of its two first derivatives along the 
real axis. Along this approximate path, the integrand does not oscillate 
too much, allowing smooth numerical integration, whose result is drawn as 
continuous lines in Fig. [2]. 

We extended the curves down to unrealistically small probabilities, in 
order to illustrate the peculiarities of the distribution, both in shape and in 
size dependence. Compared to the right hand side of the peak, the curve 
on the left is much broader and it broadens substantially as L increases. A 
similar behaviour of the Ising magnetization distribution at zero field has 
been analyzed in |L0]]. We now correlate these features to the singular part 
of <j>(z), via an analytical estimate of P / g t;J 4(e) valid at e < 0. 



The integral (|4q ) is also equal to that of minus twice the imaginary part 
of the same integrand above the cut: 

Pf3 t a(c) = -/ dx exp (A r (Re 4>(x) — e x)) sin(^4 r Im <f>(x)). (54) 

For — e not too small, the integrand is exponentially damped when |x| in- 
creases, the more so A r is large, while it vanishes at small |x| due to the 



behaviour (41) of Imcj). Thus for A r large, but fixed, and |x| small enough, 
Im (j) can be replaced by its approximation (^) and sin(^4 r 7m 4>) by its ar- 
gument. The existence of a saddle point x s follows. Neglecting Re(f>(x s ) in 
front of e x s , one finds 



/ 8 X 1 / 3 

*■ " -feffi) • (55) 

-^P0„A = ^(A r | e |r 7 / 3 exp(-(A.|e|) 2/3 ) , e<0. (56) 

Hence we have the remarkable scaling prediction that, up to a factor 
A 2 , the distribution depends on both the energy and the size through the 
single variable A r |e|. Moreover, the argument of the exponential damping 
of the probability is not extensive and its e-dependence is less than linear. 
These two properties explain the observed broadening of the distribution, 
and their very derivation emphasizes the role of a non vanishing imaginary 
part. Fig. [2] shows that the approximation (|56|) works quantitatively, the 
more so the product |-A r e| is large, which actually insures that x s and thus 
A r Im 4> are small. The behaviour (|56|) will be important for our discussion 
of metastability (see Section (5)). 
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4.4 A Numerical Simulation at q=9 



We have just seen that the scaling law (|5q ) sets in for A r |e| sufficiently 
large. Hence a numerical simulation designed to check it should accumulate 
reasonable statistics for events whose probability relative to the peak value 
is quite small. This means very extensive simulations, typically of the size 



of those of Refs.|l7|, IS]. However, Fig. [2] shows that substantial finite size 
effects already appear in q^^i^) ; Eq.(|50|), in a crossover region in between 
e = and the above mentioned scaling domain. In order to complement 
the comparison with simulations made for the cumulants (Section (3.1)), by 
a direct study of finite size effects in energy distributions, we performed a 
medium size simulation for q = 9, where the correlation length is about 1.4 
larger than at q = 10, and on lattices of sizes L = 80, 120, and 200. 

We used a Glauber type of algorithm. The site to be updated is chosen 
at random and its new spin value determined according to its Boltzmann 
weight. At P = Pt = log(4) we performed runs consisting for each L of about 
10 6 sets of L 2 random hits at the lattice sites. We measured the internal en- 
ergy density E every 10 sets, thus collecting a sample of ~ 10 5 configuration 
energies for each L value. Starting from a disordered configuration (spin 
chosen at random on each site) , the histogram obtained is proportional to 
the normalized probability density Pp tl A discussed above. For comparing the 
numerical data to the model distribution, Fig. (3), we divide each of them 
by its value at a reference energy which we choose to be Eq = —1, and as a 
function of the internal energy E we plot the quantity 




,./:') - l»*-]„M r> ,,, \ ) • ( 57 ) 

(58) 



Here we restored the physical variables using Eqs.(|47|, 48) with /2 = 12.8 as 
explained in Section (4.3). 

At each value of E, the error quoted is estimated from the fluctuations 
around the mean observed in subsamples of 10 3 to 10 4 measurements. In this 
range, the error was found to depend weakly upon the bin size, in agreement 
with an estimated autocorrelation length of order 10 2 for the average energy 
(i.e. in the peak region; it increases away of it). For the clarity of the figure, 
we suppressed the data points for which the statistics is too small to bring 
any useful information. 

We see on Fig. [3] that both the shape in energy and the dependence in 
size observed in the simulation are correctly reproduced by our construction, 
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an absolute prediction since the only parameter / was previously fixed from 
the study in section (3.1) of the cumulants at q > 10. There may be a 
tendancy for the size dependence to be too strong in the model. We do not 
consider it as significant enough to justify a retuning of / or /2, the more 
so we did not include in the errors shown that due to the fluctuations at 
the reference energy Eq. This error induces an overall uncertainty in the 
relative vertical positions of the curves. We cannot exclude either that finite 
size effects, in the sense of a residual L dependence in the free energy, are 
still present at the lowest L. 

We conclude that the adequacy of the model to predict quantitatively the 
strong size dependence observed in the disordered phase energy distribution 
below the peak, as well as its (nearly) size independent shape above the 
peak, constitutes a good evidence in favour of the analytical structure of the 
free energy inherent to the droplet picture. 



5 Beyond the transition temperature. Metasta- 
bility 

Metastability of the disordered phase refers to the possibility for the sys- 
tem to stay disordered above fit- Standard phenomenology a la Van der 
Waals, or Landau mean field theory, would state that the thermodynami- 
cal average energy density, which below 0t coincides with < E > d (/?) = 
can be continued up to a so-called spinodal value /3 sp , sepa- 
rating a metastable from an unstable region. 

In the droplet picture, as well as in field theoretic approaches, the point 
fit is a branch point, and analytic continuation above /9t is ill-defined. An 
important question, not considered in this paper, concerns the associated 
dynamics. An extensive review can be found in ||; see also |Xl]| . This field 
is still subject to active research in the context of Ising-like or liquid/ vapour 
transitions. After Langer j|, |3l|] and followers, the nucleation rate of a 
metastable state is proportional to the imaginary part of the free energy 
along the cut. Most of the recent results concern the Ising case (See for 



example |36|, [37], |38[). To our knowledge, these dynamical aspects have 
not been studied for the Potts model above q c . Here we limit ourselves 
to a discussion of which thermodynamical properties can be assigned the 
metastable states, starting from these features of energy distributions which 
we found in Section (4) as specifically due to the singular structure of the 
free energy. For the Ising case, related considerations can be found in [pLCf] . 
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5.1 Energy Distribution away from 0t by Reweighting 

On a finite lattice, the extrapolation of a given energy distribution from fit 
to any (larger or smaller) j3 is unambiguously defined by the corresponding 
reweighting of the Boltzmann factor. Hence, if Pp t a (E) * s the pure disor- 
dered phase distribution, then the (unnormalized) distribution at (3 is (we 
write it both in physical and rescaled variables ( fi7 , 48, 4^) 



D PtA (E) = exp(-A(E-EW)(p-p t ))Pp tA (E). (59) 
D(3 A {e) = exp(A r ez)P Pt , A (e). (60) 



The reweighting ( |59D is often used to smoothly interpolate, or slightly ex- 
trapolate numerical data taken at some fixed (3 point (here at (3t). The 



results of [17 1 show that it is feasible to perform a high statistics simula- 
tion at Pt without tunneling from one phase to another, so that an accurate 
numerical determination of the pure phase Pp t ,A {E) was obtained. 

We got this distribution for finite lattices from a theoretical ansatz for 
the thermodynamical limit of the free energy, and now get its continuation 
beyond fit by (pS|)- The generic shapes of the reweighted distributions as 
is varied are shown in Fig. [4] for q = 9 and L = 120. We comment them in 
the light of the results of Section (4). For e positive, Pp t ,A decreases faster 
than a Gaussian (section 4.2), and the only effect of reweighting there is 
to displace the peak position, whatever the sign of z is. For e and z both 
negative (see sections 4.3 and 4.4), the exponential growth of the Boltzmann 



factor eventually wins over the smoother behaviour (56) of Pp t ,A- Hence, 
under overcooling the distribution finally blows up for \z\ large enough, 
leading to a minimum at E = E m , visible in the upper curve on the left of 
Fig. [4], which exhibits the cubic like shape typical of a metastable situation. 
In between E m and the location Em of the maximum, all curves have an 
inflexion point at the same value E = E*: It is the inflexion point of the 
(finite size) microcanonical entropy Sa(E). 

The depth of the minimum relative to the peak height is a measure 
of the barrier to be crossed by the system in order to flip to its stable 
state, and thus controls the lifetime of the metastable state. There is no 
meaning in considering the distribution ( |59| ) below E m , where for the least 
the contributions of the ordered states, obtainable from ours by duality, and 
of mixed phase configurations [|TJ] must be included. As a definition, we 
take the distribution (|59|) cut below E m as representative of the metastable 
disordered phase. 



19 



On a given lattice, the barrier disappears when [3 becomes equal to /?*, 
the slope at E* of 

dp, A (E) = jLog(Dp >A (E)). (61) 

The size dependent point (/?*, E*) A , which plays the role of the classical 
spinodal point, is thus determined by: 

^d PA {E)\ E=Et = 0, (62) 

^dp A (E) {E=Et = /?*-/?. (63) 

At L=120 and q=9, again chosen for illustration, we find from the model 
distribution at fit (the dotted curve of Fig. [4] ) 

E* = -1.104 (64) 
P* = fy. + 3.44 10" 3 = 1.38973 . (65) 

Our statistics at /3t (Section (4.4)) are not sufficient to give significant results 
by reweighting, and we explored the above issue by new simulations in the 
vicinity of fit- For each /3 above /3t, where flips from the d to the o phase may 
become frequent, we performed ten independent runs, each run being ten 
times shorter than at (3f Flips, seen as jumps of the lattice average energy 
density El, do occur, but the time histories show that the configurations 
with El larger than ~ -1.24 can be safely labelled as disordered. Although 
an accurate determination of E* from the data is not possible , we consis- 
tently find that the distribution dg A (E) is essentially a straight line around 
E = —1.1, and we measure its slope ( |63| ) there. Our numerical results as a 
function of ft are reported in Fig. [5], and compared with the prediction (63) 



with (3* fixed by the model (p5|). The agreement is significant: the spinodal 
point location could be identified numerically as the value of (3 where data 
points for the left hand side of (|63|) linearly extrapolate to 0. 

We do not repeat the whole analysis for other lattice sizes. The impor- 
tant point is that the slope of dg t A (E) in the relevant energy region decreases 
as L increases, as seen on Fig. [3]. Our discussion of Section (4.3) interpreted 
this size dependent behaviour of the energy distribution in terms of the sin- 
gularity structure of the free energy. In particular, Eq. ( |56[ ) indicates that, 
as A — ► oo, dg >A (E) as well as its derivative eventually goes to 0, at least 
for any E held fixed sufficiently far below the peak. Via Eq. (|63|), this in 
turn implies that (3* (A) — > /9t when A — > oo. Hence we conclude that the 
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metastability interval (5* — (3t shrinks to zero in the thermodynamical limit, 
reproducing the result of the Maxwell construction applied to Van der Waals 
theory. 



5.2 Free Energy for a Finite System in a Metastable State 

We have defined the disordered energy distribution up to j3* as the dis- 
tribution ( |59| ) cut below its minimum at E m . It is then natural to use it 
in the standard relations (jH, |5|) to associate a free energy F A {(3) to the 
metastable disordered phase. So for /?</?*, we set 

J^g eMA(F A (P)-F A ((3 t ))=Y / 'D^ A (E), (66) 

4> A {z) = (y) 3 ^ ( f a(P) - F (Pt) - hip - A))- (67) 

In the above equations J2 E ' means summation above max [E m , —2]. Note 
the index A in F A and <j) A , to be explained soon. For consistency of course, 
one must find that, for z > 0, <j) A coincides with (ft, which by construction is 
A independent. This can be checked by plugging in the representation 



in the reweighted distribution (|^), and the result in (36, 67), leading to 
exp(^4 r (f) A (z)) = — — / de / dp exp [A r ((f)(z + ip) - e (z - z + ip))] . 

2 7T J em J-oo 

(68) 

The lower bound e m corresponds to the cut E m , when necessary. For z > 0, 
one may choose z = z, e m can be set to — oo with an exponentially small 
error and integration over e yields the desired result. 

The choice z = z cannot be done for z* < z < 0, where z* corresponds 
to 0* via (f49[), due to the branch point at z = 0. Furthermore, the cut at e m 
becomes relevant. For both reasons, (fi A as defined by ( |67| ) actually becomes 
A-dependent. We illustrate this definition of the free energy for metastable 
states by the continuous line of Fig. [6], drawn for L = 120 and q = 9. Also 
shown is the result of the integral (^) with e m replaced by a fixed cut off 
(here corresponding to a cut at E cu t=-1.2A). For completness, we also plot 
the real part of the input 4>{z) ( |37|) for z < (dashed line). All definitions 
provide the same result for z > 0. 

Let us add a few comments. First one observes that <j) A and Real(/> are 
quite similar just below z = 0. We know that Im</> is very small there, 
or equivalently, that the droplet expansion (|i~4D can be safely troncated 
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at the size t c where the term of order i is minimum [Q, [l|, J3f|. Near the 
end point z*, one then notice the intriguing fact that <j)A is not convex. 
However, because the restriction on the sum fl66|) depends on [3 (via E m ), the 
energy cumulants do not coincide with the successive derivatives of Fa(/3), 
they must be computed by actual averages over the distribution (^). In 
particular, the quadratic fluctuations of the energy, of course always positive, 
are not proportional to <p A , which thus may be negative. 

Likewise, the fluctuation-dissipation theorem does not hold and the spe- 
cific heat, C = —[3 2 d/df3 < E > is proportional neither to the quadratic 
fluctuations of the energy, nor to the second derivative of (ft a ■ Averaging the 
energy, we compute < E > and find it has a minimum near (3* , showing that 
C (after a maximum) does become negative. Such a phenomenon disappears 
in the thermodynamical limit since, in this limit, (3* — > (3f If a fixed cut is 
set (dots in Fig. (5)), one recovers the familiar connection between energy 
cumulants and derivatives of the free energy, and <j> is convex. Our choice of 
cutting at E m looks physically sensible. On a finite lattice, it leads to a situ- 
ation similar to that of Van der Waals theory. In the thermodynamical limit, 
it defines a disordered microcanonical entropy density S^ d \E), defined and 
concave for E^ < E < 0. A similar construction of its ordered counterpart 
S^°'(E) is immediate through duality for —2 < E < E^°\ The entropy S(E) 
for the whole system is finally achieved as usual by supplementing S^°\E) 
and SW{E) by their common tangent straight line between E^ and E^ d \ 
In this way, our approach provides a model entropy S(E) for any E. 



6 Summary and Conclusions 

We have applied a Fisher's version of the droplet picture to the 2D g-states 
Potts model in the vicinity of the temperature driven transition, which is 
first order above q c = 4. The droplet parameters were constrained by re- 
quiring that critical properties of the free energy are recovered as q — > q c , 
complying with Fisher's scaling ideas and previous empirical observations 
based on large q expansions of the energy cumulants [|^, ^0|, 21]. Using du- 



ality of the Potts model, we focussed on the disordered phase properties and 
showed that, at the price of a unique (/-independent free scale, the energy 
cumulants of order <10 measured in a numerical simulation by Janke and 
Kappler were accurately reproduced (Fig.[l]). The analytical structure of 
the thermodynamical free energy in the inverse temperature (3 plane was 
then analyzed, and a branch point essential singularity at the transition 
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fully characterized. Next we investigated the consequences of these prop- 
erties directly on the energy distribution currently generated in numerical 
simulations on finite lattices. In particular, evaluating this distribution an- 
alytically far away from its peak, we established a direct link between its 
size dependence there and the free energy discontinuity around the branch 
cut. Finite size effects are also present close to the peak of the distribution 
and in good agreement with the result of a numerical simulation which we 
performed at q = 9 for comparison. We finally discussed static aspects of 
metastability. We defined a free energy beyond the transition point fit via 
the reweighting of the energy distribution at finite size, showing that, again 
due to the essential singularity, this can be done up to a size dependent spin- 
odal point only, which moreover coincides with fit in the thermodynamical 
limit. 

We noted earlier that various versions of the generic droplet picture 
have been used to study the first order transition of the Ising model below 
its critical temperature, and the associated issues of singularity of the free 
energy, metastability, and nucleation rates. As already mentioned in Section 
(2.3), the Fisher's version which we apply as been questionned for the Ising 
case in that it assigns a droplet an interfacial energy which grows faster 
than the perimeter (cr > 1/2 in (14)), which sounds odd at least for large 
droplets. Binder has shown |5], [| how scaling prescriptions valid close to the 
critical point (T = T c , h = for Ising) can reconcile the Fisher's description 
and the original geometric picture of nucleation. We expect that similar 
arguments apply to the Potts case around q = q c , (3 = fit- In any case, our 
work shows that following the Fisher's point of view provides an extremely 
efficient and economical parametrization of the pure phase free energies, 
as particularly illustrated by its adequacy to describe the ten first energy 
cumulants: a) their measurements in numerical simulations |L7]] require the 
predicted, non-geometrical, values a = 2/3 andr = 7/3 (see Eqs.(31, 32)), b) 
the pseudo interfacial tension u scales as prescribed by the Widom exponent 
for the conventional tension. 

We see three directions in which it would be interesting to pursue with 
the Potts model, noting that the above indices differ from their geometrical 
values more than they do for the Ising case, i) Derivation of exact results, 
in analogy with what has been done for Ising or Ising-like models on clus- 
ter distributions, large order cumulants, metastable states, ii) Dynamics at 
and near the transition and relationship between nucleation rates and sin- 
gularities of the free energies, iii) Relationship between finite size properties 
of energy distributions at the transition and the discontinuity of the pure 
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phase free energy. The existence of such a connection is demonstrated by our 
derivation of Eq. ( |56|) (which can be easily extended to arbitrary values of the 
indices a and r). Accordingly, high statistics numerical data |17], [l8| may 
provide a direct access to the free energy discontinuities, without requiring 
heavy dynamical investigations. 
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Figure 1: Comparison of the prediction of Eq. (^8|) (continuous line) with the 
numerical data of Ref.[]l7|], versus the order n. The large squared symbols 
at n = 9 and 10 indicate orders of magnitude only. The scale factor S n , 
which contains the only free parameter /, exhausts the full q-dependence of 
the / n 's. 
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Figure 2: The energy distribution of the disordered phase at q = 9 as 
a function of e for L = 80, 120, 200. The continuous curves result from 
numerical integration of ( [46|) and the symbols represent analytical approx- 
imations (see text). As L increases, the strong broadening for e < is a 
manifestation of the free energy singularity at j3t- 
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Figure 3: Comparison, for pp u A(E), Eq.(|57|), between the model predictions 
and numerical data taken at q = 9, L=80, 120, and 200 (from bottom to 
top). 
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Figure 4: For the disordered phase at q = 9, L = 120, energy distribu- 
tions for ((3 — fit) varying from —2 10~ 3 to 2 10~ 3 . On the left of the figure, 
temperature decreases from the bottom to the top curves (overcooling) , the 
dotted curve corresponding to f3 = (3f The crosses correspond to the inflex- 
ion points which all occur at the same E* = —1.104. As the temperature 
is lowered below the transition, a minimum of the distribution enters the 
graph from the left. 
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Figure 5: The slope A &&,A (E)\ E - E * as a function of (3 from numerical 
simulations (o). The straight line represents Eq. (|63|), with (3* given by the 
model fl65|): the slope vanishes at the spinodal value (3* (□). 



31 



0.8 



0.6 



0.4 



0.2 



0.0 



^ i i r 



l i i r 



L=120 



q=9 




i i r 



1.0 



0.5 



0.0 



0.5 



Figure 6: The rescaled free energy for metastable states (z < 0). Are shown 
the sum ( |67[ ) where the cut e m is the z dependent minimum of the reweighted 
distribution (continuous curve), the same sum with a fixed cut (E cu t=-1.24, 
dots), and the real part of the input (j)(z) (dashed curve). The symbol □ 
marks the end point of the metastability region. It moves right and down 
with increasing L, ending at the origin in the thermodynamical limit. All 
definitions give the same answer (j)(z) for z > 0. 
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